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Abstract 

Let H 0 , ..., H n be m x in matrices with entries in Q and Hankel structure, i.e. 
constant skew diagonals. We consider the linear Hankel matrix H(x.) = Ho + x\Hi + 
■ ■ ■ + x n H n and the problem of computing sample points in each connected compo¬ 
nent of the real algebraic set defined by the rank constraint rank(H(x)) < r, for a 
given integer r < m — 1. Computing sample points in real algebraic sets defined by 
rank defects in linear matrices is a general problem that finds applications in many 
areas such as control theory, computational geometry, optimization, etc. Moreover, 
Hankel matrices appear in many areas of engineering sciences. Also, since Hankel 
matrices are symmetric, any algorithmic development for this problem can be seen 
as a first step towards a dedicated exact algorithm for solving semi-definite program¬ 
ming problems, i.e. linear matrix inequalities. Under some genericity assumptions 
on the input (such as smoothness of an incidence variety), we design a probabilistic 
algorithm for tackling this problem. It is an adaptation of the so-called critical point 
method that takes advantage of the special structure of the problem. Its complexity 
reflects this: it is essentially quadratic in specific degree bounds on an incidence 
variety. We report on practical experiments and analyze how the algorithm takes 
advantage of this special structure. A first implementation outperforms existing im¬ 
plementations for computing sample points in general real algebraic sets: it tackles 
examples that are out of reach of the state-of-the-art. 
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1 Introduction 


Problem statement and motivation Let Q, R, C be respectively the fields of rational, 
real and complex numbers, and let m, n be positive integers. Given m x m matrices 
Ho, Hi,..., H n with entries in Q and Hankel structure, i.e. constant skew diagonals, we 
consider the linear Hankel matrix R(x) = H 0 + X\H\ + ... + x n H n , denoted H for short, 
and the algebraic set 

H r = {x 6 C" : rankR(x) < r}. 

The goal of this paper is to provide an efficient algorithm for computing at least one 
sample point per connected component of the real algebraic set H r C R n . 

Such an algorithm can be used to solve the matrix rank minimization problem for H. 
Matrix rank minimization mostly consists of minimizing the rank of a given matrix whose 
entries are subject to constraints defining a convex set. These problems arise in many 
engineering or statistical modeling applications and have recently received a lot of atten¬ 
tion. Considering Hankel structures is relevant since it arises in many applications (e.g. 
for model reduction in linear dynamical systems described by Markov parameters, see [16, 
Section 1.3]). 

Moreover, an algorithm for computing sample points in each connected component of 
H r fl R n can also be used to decide the emptiness of the feasibility set S = {x 6 R n : 
H(x) >z 0}. Indeed, considering the minimum rank r attained in the boundary of S, it 
is easy to prove that one of the connected components of H r fl R™ is actually contained 
in S. Note also that such feasibility sets, also called Hankel spectrahedra, have recently 
attracted some attention (see e.g. [3]). 

The intrinsic algebraic nature of our problem makes relevant the design of exact algorithms 
to achieve reliability. On the one hand, we aim at exploiting algorithmically the special 
Hankel structure to gain efficiency. On the other hand, the design of a special algorithm 
for the case of linear Hankel matrices can bring the foundations of a general approach to 
e.g. the symmmetric case which is important for semi-definite programming, i.e. solving 
linear matrix inequalities. 


Related works and state-of-the-art Our problem consists of computing sample 
points in real algebraic sets. The first algorithm for this problem is due to Tarski but its 
complexity was not elementary recursive [23]. Next, Collins designed the Cylindrical Al¬ 
gebraic Decomposition algorithm [5]. Its complexity is doubly exponential in the number 
of variables which is far from being optimal since the number of connected components 
of a real algebraic set defined by n-variate polynomial equations of degree < d is upper 
bounded by 0(d) n . Next, Grigoriev and Vorobjov [12] introduced the first algorithm 
based on critical point computations computing sample point in real algebraic sets within 
d°( n ) arithmetic operations. This work has next been improved and generalized (see [2] 
and references therein) from the complexity viewpoint. We may apply these algorithms 
to our problem by computing all (r + l)-minors of the Hankel matrix and compute sample 
points in the real algebraic set defined by the vanishing of these minors. This is done in 
time ((^jJ + r°( n ) however since the constant in the exponent is rather high, 

these algorithms did not lead to efficient implementations in practice. Hence, another se- 
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ries of works, still using the critical point method but aiming at designing algorithms that 
combine asymptotically optimal complexity and practical efficiency has been developed 
(see e.g. [1, 19, 11] and references therein). 

Under regularity assumptions, these yield probabilistic algorithms running in time which 
is essentially 0(d 3n ) in the smooth case and 0(d 4n ) in the singular case (see [18]). Practi¬ 
cally, these algorithms are implemented in the library RAGlib which uses Grobner bases 
computations (see [8, 22] about the complexity of computing critical points with Grobner 
bases). 

Observe that determinantal varieties such as l~L r are generically singular (see [4]). Also 
the aforementioned algorithms do not exploit the structure of the problem. In [14], we 
introduced an algorithm for computing real points at which a generic linear square matrix 
of size m has rank < m — 1, by exploiting the structure of the problem. However, because 
of the requested genericity of the input linear matrix, we cannot use it for linear Hankel 
matrices. Also, it does not allow to get sample points for a given, smaller rank deficiency. 


Methodology and main results Our main result is an algorithm that computes sam¬ 
ple points in each connected component of "H r fl M n under some genericity assumptions on 
the entries of the linear Hankel matrix H (these genericity assumptions are made explicit 
below). Our algorithm exploits the Hankel structure of the problem. Essentially, its com¬ 
plexity is quadratic in a multilinear Bezout bound on the number of complex solutions. 
Moreover, we find that, heuristically, this bound is less than j (”^ r ) . Hence, for 

subfamilies of the real root finding problem on linear Hankel matrices where the maximum 
rank allowed r is fixed, the complexity is essentially in {nm)°^ r \ 

The very basic idea is to study the algebraic set T-L r C C n as the Zariski closure of the 
projection of an incidence variety, lying in C n+r+1 . This variety encodes the fact that the 
kernel of H has dimension >m — r. This lifted variety turns out to be generically smooth 
and equidimensional and defined by quadratic polynomials with multilinear structure. 
When these regularity properties are satisfied, we prove that computing one point per 
connected component of the incidence variety is sufficient to solve the same problem for 
the variety di r 0 M”. We also prove that these properties are generically satisfied. We 
remark that this method is similar to the one used in [14], but in this case it takes strong 
advantage of the Hankel structure of the linear matrix, as detailed in Section 2. This also 
reflects on the complexity of the algorithm and on practical performances. 

Let C be a connected component of "H r flM”, and and Hi, 7Ti be the canonical projections 
n x : (aq,..., x n , y 1: ... , y r +i) —> aq and 7Ti : (aq,..., x n ) —> aq. We prove that in generic 
coordinates, either (i) ^(C) = R or (ii) there exists a critical point of the restriction of 
n x to the considered incidence variety. Hence, after a generic linear change of variables, 
the algorithm consists of two main steps: (i) compute the critical points of the restriction 
of n x to the incidence variety and (ii) instantiating the first variable aq to a generic value 
and perform a recursive call following a geometric pattern introduced in [19]. 

This latter step (i) is actually performed by building the Lagrange system associated to 
the optimization problem whose solutions are the critical points of the restriction of 7iq to 
the incidence variety. Hence, we use the algorithm in [15] to solve it. One also observes 
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heuristically that these Lagrange systems are typically zero-dimensional. 

However, we were not able to prove this finiteness property, but we prove that it holds 
when we restrict the optimization step to the set of points x G l~L r such that rank H (x) = p, 
for any 0 < p < r. However, this is sufficient to conclude that there are finitely many 
critical points of the restriction of to l~L r fl M n , and that the algorithm returns the 
output correctly. 

When the Lagrange system has dimension 0, the complexity of solving its equations is 
essentially quadratic in the number of its complex solutions. As previously announced, by 
the structure of these systems one can deduce multilinear Bezout bounds on the number of 
solutions that are polynomial in nm when r is fixed, and polynomial in n when m is fixed. 
This complexity result outperforms the state-of-the-art algorithms. We finally remark that 
the complexity gain is reflected also in the first implementation of the algorithm, which 
allows to solve instances of our problem that are out of reach of the general algorithms 
implemented in RAGlib. 


Structure of the paper The paper is structured as follows. Section 2 contains prelim¬ 
inaries about Hankel matrices and the basic notation of the paper; we also prove that our 
regularity assumptions are generic. In Section 3 we describe the algorithm and prove its 
correctness. This is done by using preliminary results proved in Sections 5 and 6. Section 
4 contains the complexity analysis and bounds for the number of complex solutions of 
the output of the algorithm. Finally, Section 7 presents the results of our experiments on 
generic linear Hankel matrices, and comparisons with the state-of-the-art algorithms for 
the real root finding problem. 


2 Notation and preliminaries 

Basic notations We denote by GL(n, Q) (resp. GL(n, C)) the set ofnxn non-singular 
matrices with rational (resp. complex) entries. For a matrix M G C mxm and an integer 
p < m, one denotes with minors (p, M ) the list of determinants of p x p sub-matrices of 
M. We denote by M' the transpose matrix of M. 

Let Q[x] be the ring of polynomials on n variables x = (aq,. .., x n ) and let f = (/i,..., f p ) G 
Q[x] p be a polynomial system. The common zero locus of the entries of f is denoted by 
Z(f) C C n , and its dimension with dim Z(f). The ideal generated by f is denoted by (f), 
while if V C C n is any set, the ideal of polynomials vanishing on V is denoted by /(V), 
while the set of regular (resp. singular) points of V is denoted by reg V (resp. sing V). 
If f = (fi ,..., f p ) C Q[x], we denote by Df = ( dfi/dxj ) the Jacobian matrix of f. We 
denote by reg(f) C Z(f) the subset where D f has maximal rank. 

A set £ C C n is locally closed if £ = Z fl & where Z is a Zariski closed set and O is a 
Zariski open set. 

Let V = Z( f) C C" be a smooth equidimensional algebraic set, of dimension d, and let 
g: C n —> C p be an algebraic map. The set of critical points of the restriction of g to 
V is the solution set of f and of the (n — d + p)— minors of the matrix D (f, g), and it 
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is denoted by crit (g, V). Finally, if £ C V is a locally closed subset of V, we denote by 
crit (g, £) = £ fl crit (g, V). 

Finally, for M £ GL(n,C) and / £ Q[x], we denote by / M (x) = /(Mx), and if f = 
(fi, ■ ■ ■, / P ) C Q[x] and V = Z( f), by V M = Z(f M ) where i M = (/ X M , ..., / p M ). 

Hankel structure Let {hi,..., hzm-i} C Q. The matrix H = G Q mxm 

is called a Hankel matrix, and we use the notation H = Hankel (hi,...,/i 2 m-i)- The 
structure of a Hankel matrix induces structure on its kernel. By [13, Theorem 5.1], one 
has that if H is a Hankel matrix of rank at most r, then there exists a non-zero vector 
y = (yi ,..., y r +i ) £ Q r+1 such that the columns of the m x (m — r ) matrix 

> 0 ... O' 

r(y)= 0 y : 

: ••• 0 

.0 ... 0 y 

generate a (m — r)— dimensional subspace of the kernel of H. We observe that HY(y) is 
also a Hankel matrix. 

The product HY(y) can be re-written as a matrix-vector product Hy , with H a given 
rectangular Hankel matrix. Indeed, let H = Hankel(hi,..., h^m-x)■ Then, as previously 
observed, HY(y) is a rectangular Hankel matrix, of size m x (m — r), whose entries 
coincide with the entries of 


Hy = 


hi ■ ■ ■ h r +1 


2/i 

^2m-r-l • • • h'lrn— 1 


Vr+ 1_ 


Let H(x) be a linear Hankel matrix. From [6, Corollary 2.2] we deduce that, for p < r, 
then the ideals (minors (p + 1, H(x))) and ^minors (p + 1, H(x)) S J coincides. One deduces 
that x = (# 1 ,..., x n ) £ C n satisfies rank H(x) — p if and only if it satisfies rank iL(x) = p. 


Basic sets We first recall that the linear matrix H(x) = Hq + X\Hi + ... + x n H n , where 
each Hi is a Hankel matrix, is also a Hankel matrix. It is identified by the (2 m — l)(n+ 1) 
entries of the matrices Hi. Hence we often consider H as an element of C ( ^ 2m_1 ^ n + 1 ). p or 
M £ GL(n, Q), we denote by iL M (x) the linear matrix H(Mx). 

We define in the following the main algebraic sets appearing during the execution of our 
algorithm, given H £ 0 < p < r, M £ GL(n, C) and u = {u\,... ,u p+ 1 ) £ 

Qp +1 . 

Incidence varieties. We consider the polynomial system 

f(tf M ,u,p) : C n x C p+1 —> C 2 ™-?- 1 x C 

(x,y) ^ [(H(M x) y)', u'y — l) 
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where H has been defined in the previous section. We denote by J(H M , u, p) = Z(f p (H M , u)) C 
C" +p+ i and simply J = JfiH M , u ,p) and f = f (H M , u ,p) when p, H, Ad and u are clear. 

We also denote by JC(H M , u ,p) = J(H M , u ,p) D {(x, y) G C n+P+1 : rank //(x) = p}. 

Fibers. Let a G Q. We denote by u,p) (or simply f Q ) the polynomial system 

obtained by adding X\ — a to f (H M , u ,p). The resulting algebraic set Z(f a ), denoted by 
J a , equals J fl Z{x\ — a). 


Lagrange systems. Let v G Q 2m ~ p . Let D \f denote the matrix of size c x (n + p) obtained 
by removing the first column of D f (the derivative w.r.t. aq), and define £ = £(H M , u, v, p) 
as the map 


^n+2m+l 

(x,y,z) 


_ y (Qr).+2»n+l 

^ (H(Ix)y,u'y 


1, z'D if, v'z 


1 ) 


where z = ..., Z 2 m - P ) stand for Lagrange multipliers. We finally define Z(H M , u, v, p ) = 

Z(£(H M ,u,v,p)) C C n+2m+1 . 


Regularity property G We say that a polynomial system f G Q[a;] c satisfies Property 
G if the Jacobian matrix D f has maximal rank at any point of Z( f). We remark that 
this implies that: 

1. the ideal /(f) is radical; 

2. the set Z( f) is either empty or smooth and equidimensional of co-dimension c. 

We say that £(H M , u, v, p) satisfies G over /C(// M , u, p) if the following holds: for (x, y, z) G 
Z(H m , u, v,p) such that (x, y) G JC(H M , u ,p), the matrix D (/(ZZ M , u, v, p)) has maximal 
rank at (x, y, z). 

Let u G Q p+1 . We say that H G c( 2m_1 k n+1 ) satisfies Property G if f(H,u,p) satisfies 
Property G for all 0 < p < r. 

The first result essentially shows that G holds for f(// M ,u,p) (resp. f a (H M , u,p)) when 
the input parameter H (resp. a) is generic enough. 

Proposition 1 Let Ad G GL(n, C). 

(a) There exists a non-empty Zariski-open set C C^-Tin+i) suc p that, if H G 
Jdf fl Q( 2m_1 )( n + 1 ) ? f or all 0 < p < r and u G Q p+1 — {0}, i(H M ,u,p ) satisfies 
Property G; 

(b) for H G iff, and 0 < p < r, if J{H M , u ,p) 0 then dim TL P < n — 2m + 2p+l; 

(c) For 0 < p < r and u G Q p+1 , if i(H M ,u,p) satisfies G. there exists a non-empty 
Zariski open set sZ cC such that, if a G stf, the polynomial system f Q satisfies G; 

Proof : Without loss of generality, we can assume that M = I„. We let 0 < p < r, 
u G Q p+1 — {0} and recall that we identify the space of linear Hankel matrices with 
C (2m-I)(n+1)_ This space is endowed by the variables t)k,e with 1 < k < 2m — 1 and 
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0 < £ < n; the generic linear Hankel matrix is then given by = Sj 0 + X\$)i + • • • + x n S) n 
with fii = Hankel(t)i j,..., 

We consider the map 

q ■ £n+(p+l)+(2m-l)(n+l) _ y (£2m—p 

(x,y ,H) i—> f(H,n,p) 

and, for a given H G C^ 2m_1 ^ n+1 \ its section-map Qh '■ C n+ ^ p+1 ) —» C 2m_p sending (x, y) 
to g(x, y, H). We also consider the map q which associates to (x, y, H) the entries of H y 
and its section map qu\ we will consider these latter maps over the open set O — {(x, y) G 
C n+P+1 | y 7 ^ 0}. We prove below that 0 is a regular value for both q H and q H . 

Suppose hrst that g _ 1 (0) = 0 (resp. g _ 1 (0)). We deduce that for all H G (CJ 2m_1 ^ n+1 ), 
9 ^( 0 ) = 0 (resp 5 ^( 0 ) = 0) and 0 is a regular value for both maps Qh and qn- Note 
also that taking Jrf? = C( 2m_1 )( n+1 ), we deduce that f(i7, u,p) satisfies G. 

Now, suppose that g _ 1 (0) is not empty and let (x, y , H) G g _ 1 (0). Consider the Jacobian 
matrix D q of the map q with respect to the variables x, y and the entries of JJ, evaluated 
at (x, y, H). We consider the submatrix of D q by selecting the column corresponding to: 

• the partial derivatives with respect to f)i >0 ,..., 

• the partial derivatives with respect to yi ,..., y p+ 1 . 

We obtain a (2m — p) x (2m + p) submatrix of D q: we prove below that it has full rank 
2 771 — p. 

Indeed, remark that the 2m — p — 1 first lines correspond to the entries of Hy and last 
line corresponds to the derivatives of u'y — 1. Hence, the structure of this submatrix is 
as below 

~yi ... y p+ 1 0 ... 0 0 • • • O' 

0 y 1 ... y p+1 ... 0 

: yi ■■■ y P +1 0 0 

_ 0 • • • • • • 0 Ui ■ ■ ■ u p+ 1 _ 

Since this matrix is evaluated at the solution set of u'y — 1 = 0, we deduce straightfor¬ 
wardly that one entry of u and one entry of y are non-zero and that the above matrix is 
full rank and that 0 is a regular value of the map q. 

We can do the same for D q except the fact that we do not consider the partial derivatives 
with respect to 7 / 1 ,..., 7/ p+ i. The (2m— p— 1) x (2m — 1) submatrix we obtain corresponds 
to the upper left block containing the entries of y. Since q is defined over the open set O 
in which y ^ 0, we also deduce that this submatrix has full rank 2m — p — 1. 

By Thom’s Weak Transversality Theorem one deduces that there exists a non-empty 
Zariski open set M > p C C^ 2m_1 ^ n+1 ^ such that if H G J^ p , then 0 is a regular value of 
qn (resp. qu). We deduce that for H G Jrf? P) the polynomial system i(H,u,p) satisfies 
G and using the Jacobian criterion [7, Theorem 16.19], J(H, u, p) is either empty or 
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smooth equidimensional of dimension n — 2m + 2p + 1. This proves assertion (a), with 
J^ = f] 0 < p <r^p- 

Similarly, we deduce that qJj X (0) is either empty or smooth and equidimensional of di¬ 
mension n — 2m + 2p + 2. Let fl x be the canonical projection (x, y) —> x; note that for 
any x G the dimension of fl x 1 (x) fl qJj X (0) is > 1 (by homogeneity of the y-variables). 
By the Theorem on the Dimension of Fibers [21, Sect.6.3,Theorem 7], we deduce that 
n — 2m + 2p + 2 — dim('Hp) > 1. We deduce that for H G Jtf, dim (H p ) <n — 2m + 2p + 1 
which proves assertion (b). 

It remains to prove assertion (c). We assume that i(H,u,p) satisfies G. Consider the 
restriction of the map n x : C n+P+1 —» C, II 1 (x, y) = aq, to J(H,x\,p), which is smooth 
and equidimensional by assertion (a). 

By Sard’s Lemma [20, Section 4.2], the set of critical values of the restriction of Lb to 
J{H, u,p) is finite. Hence, its complement sZ C C is a non-empty Zariski open set. We 
deduce that for a G &/, the Jacobian matrix of f a (H, u ,p) satisfies G. □ 


3 Algorithm and correctness 

In this section we present the algorithm, which is called LowRankHankel, and prove its 
correctness. 

3.1 Description 

Data representation The algorithm takes as input a couple ( H , r), where H = (Hq, Hi, , H n ) 
encodes m x m Hankel matrices with entries in Q, defining the linear matrix i7(x), and 
0 < r < m — 1. 

The output is represented by a rational parametrization, that is a polynomial system 

q = (qo(t), qi{t), • • •, q n {t), q(t )) C Q [t] 
of univariate polynomials, with gcd(q,qo) = 1. The set of solutions of 

%i ~ qi{t)/q 0 {t) = 0, i = 1... 7i q(t) = 0 

is clearly finite and expected to contain at least one point per connected component of 
the algebraic set H r D M”. 


Main subroutines and formal description We start by describing the main subrou¬ 
tines we use. 

ZeroDimSolve. It takes as input a polynomial system defining an algebraic set Z C 
C n+fc and a subset of variables x = (aq,..., x n ). If Z is finite, it returns a rational 
parametrization of the projection of Z on the x-space else it returns an empty list. 

ZeroDimSolveMaxRank. It takes as input a polynomial system f = (/j,..., f c ) such that 
Z = {x G C n+k | rank(ZJf(x)) = c} is finite and a subset of variables x = (aq,... ,x n ) 
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that endows C n . It returns fail: the assumptions are not satisfied if assumptions are not 
satisfied, else it returns a rational parametrization of the projection of Z on the x-space. 

Lift. It takes as input a rational parametrization of a finite set Z C and a number 
a G C, and it returns a rational parametrization of {(a, x) : x G Z}. 

Union. It takes as input two rational parametrizations encoding finite sets Z \, Z 2 and it 
returns a rational parametrization of Z\ U Z 2 . 

ChangeVariables. It takes as input a rational parametrization of a finite set Z C C N and 
a non-singular matrix M G GL(7V, C). It returns a rational parametrization of Z AI . 

The algorithm LowRankHankel is recursive, and it assumes that its input H satisfies Prop¬ 
erty G. 

LowRankHankel(i/, r): 

1. If n < 2 m — 2r — 1 then return []. 

2. Choose randomly M G GL(n, Q), a G Q and u p G Q p+1 , v p G Q 2m ~ p for 0 < p < r. 

3. If n — 2m — 2r — 1 then return ZeroDimSolve(f(iJ, u r , r), x)). 

4. Let P = ZeroDimSolve(T(f (if, u r , r), v)) 

5. If P = [] then for p from 0 to r do 

(a) P' = ZeroDimSolveMaxRank(f(if M , u p , v p ), x); 

(b) P = Union(P, P') 

6. Q = Lift(LowRankHankel(Subs(o:i = a, H M ),r),a); 

7. return(ChangeVariables(Union(Q, P), 

3.2 Correctness 

The correctness proof is based on the two following results that are proved in Sections 5 
and 6. 

The first result states that when the input matrix H satisfies G and that, for a generic 
choice of M and v, and for all 0 < p < r, the set of solutions (x, y, z) to £(H M , u, v,p) at 
which rankiL(rr) — p is finite and contains crit (ni, IC(H m , u,p)). 

Proposition 2 Let be the set defined in Proposition 1 and let H G M J and u G 
Qp+i _ j or o < p < r. There exist non-empty Zariski open sets C GL(n, C) and 
"V C C 2rn ~ p such that if M G .M\ D Q nxn and v G Y D Q 2m ~ p , the following holds: 

(a) £(H M ,u,v,p) satisfies G over IC(H M ,u,p); 

(b) the projection of reg (£(H M , u, v,p)) on the (x, y )-space contains crit (Ilx, IC(H AI , u,p)) 


9 


Proposition 3 Let H G Jrff , 0 < p < r and d p = n — 2m + 2p + 1 and C be a connected 
component o/'H p nR n . Then there exist non-empty Zariski open sets ^ 2 C GL(n, C) and 
c C p+1 such that for any M G ^ 2 H Q nxn , ugfn Q p+1 , the following holds: 

(a) for i — 1,..., d p , 7r i{C M ) is closed; 

(b) for any a G R in the boundary ofni(C M ), 7t 1 “ 1 (q!) D C m is finite; 

(c) for any x G 7rf 1 (o!) D C M and p such that rank// p (x) = p, there exists (x, y) G 
R” x R p+1 such that (x, y) G J(H M , u,p) . 

Our algorithm is probabilistic and its correctness depends on the validity of the choices 
that are made at Step 2. We make this assumption that we formalize below. 

We need to distinguish the choices of M, u and v that are made in the different calls of 
LowRankHankel; each of these parameter must lie in a non-empty Zariski open set defined 
in Propositions 1, 2 and 3. 

We assume that the input matrix H satisfies G; we denote it by H^\ where the super 
script indicates that no recursive call has been made on this input; similarly a® denotes 
the choice of a made at Step 2 on input H^°\ Next, we denote by the input of 
LowRankHankel at the i-Wi recursive call and by s/® C C the non-empty Zariski open 
set defined in Proposition 1 applied to Note that if G we can deduce that 
satisfies G. 

Now, we denote by and the open sets defined in Propositions 1, 

2 and 3 applied to H^\ for 0 < p < r and where i is the depth of the recursion. 

Finally, we denote by g GL(n,Q), uW g Q p+1 and vj*\ for 0 < p < r, respectively 
the matrix and the vectors chosen at Step 2 of the i-th call of LowRankHankel. 

Assumption A. We say that A is satisfied if Afb), oft\ uW and satisfy: 

• g {Jt? n JZ?) nQ !Xi ; 

• a® G 

• uG G fi Q p+1 — {0}, for 0 < p < r; 

• vf) G n Q 2m_p - {0} for 0 < p < r; 

Theorem 4 Let H satisfy G. Then, if A is satisfied, LowRankHankel with input (H, r), 
returns a rational parametrization that encodes a finite algebraic set in Ti r meeting each 
connected component ofTL r flR". 

Proof : The proof is by decreasing induction on the depth of the recursion. 

When n < 2m — 2r — 1, TL r is empty since the input H satisfies G (since A is satisfied). 
In this case, the output defines the empty set. 
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When n = 2m—2r — l, since A is satisfied, by Proposition ??, either l~L r = 0 or dim T-i r = 0. 
Suppose % r = 0. Hence J r = 0, since the projection of J r on the x—space is included 
in 3i r . Suppose now that dim 7~L r = 0: Proposition 3 guarantees that the output of the 
algorithm defines a finite set containing 'H r - 

Now, we assume that n > 2 m — 2r — 1; our induction assumption is that for any i > 1 
LowRankHankel(//b\ r ) returns a rational parametrization that encodes a finite set of 
points in the algebraic set defined by rank(ifW) < r and that meets every connected 
component of its real trace. 

Let C be a connected component of "H r D M n . To keep notations simple, we denote by 
M G GL(n,Q), u p and v p the matrix and vectors chosen at Step 2 for 0 < p < r. Since 
A holds one can apply Proposition 3. We deduce that the image is closed. Then, 

either 7 t 1 (C' m ) = R or it is a closed interval. 

Suppose first that n 1 (C M ) = M. Then for a G Q chosen at Step 2, 7rf 1 (a) n C AI ^ 
0. Remark that 7rf 1 (o;) D C M is the union of some connected components of n 
R n_1 = {x = (x 2 ,...,x n ) G M n_1 : rank if ^(x) < r}. Since A holds, assertion (c) of 
Proposition 1 implies that ffd) satisfies G. We deduce by the induction assumption that 
the parametrization returned by Step 6 where LowRankHankel is called recursively defines 
a finite set of points that is contained in l~L r and that meets C. 

Suppose now that (C M ) ^ M. By Proposition 3, (C M ) is closed. Since C M is 
connected, is a connected interval, and since 7i ^ M there exists j3 in 

the boundary of ni(C M ) such that 7Ti (C m ) C [/?, +oo) or C (—oo,/3]. Suppose 

without loss of generality that 7 Ti(C' m ) C [/3, +c»), so that f5 is the minimum value attained 
by 7Ti on C M . 

Let x = ((3,x 2 ,...,x n ) G C M , and suppose that rank(if(x)) = p. By Proposition 3 
(assertion (c)), there exists y G C p+1 such that (x, y) G J(H, u,p). Note that since 
rank(if(x)) = p, we also deduce that (x, y) G JC(H, u ,p). 

We claim that there exists z G (C 2m ~ p such that (x, y, z) lies on reg (£(if M , u, v,p)). 

Since A holds, Proposition 2 implies that i(H AI , u, v,p) satisfies G over IC(H AI , u,p). Also, 
note that the Jacobian criterion implies that reg(£(H M ,u,v,p)) has dimension at most 
0. We conclude that the point x G C AI lies on the finite set encoded by the rational 
parametrization P obtained at Step 5 of LowRankHankel and we are done. 

It remains to prove our claim, i.e. there exists z G C 2m-p such that (x, y, z) lies on 
reg (£(H M , u, v,p)). 

Let C be the connected component of J(H, u,p) AI n containing (x, y). We 

first prove that /3 = ^(x, y) lies on the boundary of 7T 1 (C / ). Indeed, suppose that there 
exists (x, y) G C’ such that ^(x, y) < (3. Since C is connected, there exists a continuous 
semi-algebraic map r: [0,1] —> C with r(0) = (x, y) and r(l) = (x, y). Let p : (x, y) —> x 
be the canonical projection on the x-space. 

Note that ip o r is also continuous and semi-algebraic (it is the composition of continuous 
semi-algebraic maps), with (p o r)(0) = x, (p o r)(l) = x. Since (p o r)(6) G 7~L P for all 
6 G [0,1], then xGC. Since 7Ti(x) = 7Ti(x, y) < a we obtain a contradiction. So 7Ti(x, y) 
lies on the boundary of iri(C'). 
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By the Implicit Function Theorem, and the fact that f(H,u,p) satisfies Property G, one 
deduces that (x, y) is a critical point of the restriction of 1^! : (art,, x n , y 1, • • •, y r + 1) —> 
xi to J(H, u ,p). 

Since rank(if M (x)) — p by construction, we deduce that (x, y) is a critical point of the 
restriction of hi! to 1C(H M , u ,p) and that, by Proposition 2, there exists z G C 2m_p such 
that (x, y, z) belongs to the set reg (£(H M , u, v,p)), as claimed. □ 


4 Degree bounds and complexity 

We first remark that the complexity of subroutines Union, Lift and ChangeVariables (see [20, 
Chap. 10]) are negligible with respect to the complexity of ZeroDimSolveMaxRank. Hence, 
the complexity of LowRankHankel is at most n times the complexity of ZeroDimSolveMaxRank, 
which is computed below. 

Let (H, r) be the input, and let 0 < p < r. We estimate the complexity of ZeroDimSolveMaxRank 
with input (H M , u p , v p ). It depends on the algorithm used to solve zero-dimensional poly¬ 
nomial systems. We choose the one of [15] that can be seen as a symbolic homotopy 
taking into account the sparsity structure of the system to solve. More precisely, let 
pCQ[ii,..., x n ] and s G Q[aq ,..., x n ] such that the common complex solutions of poly¬ 
nomials in p at which s does not vanish is finite. The algorithm in [15] builds a system 
q that has the same monomial structure as p has and defines a finite algebraic set. Next, 
the homotopy system t = tp + (1 — £)q where t is a new variable is built. The system 
t defines a 1-dimensional constructible set over the open set defined by s ^ 0 and for 
generic values of t. Abusing notation, we denote by Z( t) the curve defined as the Zariski 
closure of this constructible set. 

Starting from the solutions of q which are encoded with a rational parametrization, the 
algorithm builds a rational parametrization for the solutions of p which do not cancel 
s. Following [15], the algorithm runs in time 0~(Ln°^ 88') where L is the complexity of 
evaluating the input, 8 is a bound on the number of isolated solutions of p and 8' is a 
bound on the degree of Z(t) defined by t. 

Below, we estimate these degrees when the input is a Lagrange system as the ones we 
consider. 

Degree bounds. We let ((By)', u p 'y - 1), with y = (y 1} ..., y p+1 )', defining J P (H, u p ). 

Since y ^ 0, one can eliminate w.l.o.g. y p+ 1 , and the linear form Up'y — 1, obtaining a 
system f G Q[x, y] 2m_p ~ 1 . We recall that if x (1 \ ... ,x^ c ^ are c groups of variables, and 
/ G Q[x (1 \ ..., x( c )], we say that the multidegree of / is (di ,..., d c ) if its degree with 
respect to the group x ( b i s d 3 , for j — 1 ,..., c. 

Let l = (f, g, h) be the corresponding Lagrange system, where 

(g, h) = (gi, • • •, g n - 1 , hi ,..., hp) = ZD J 

with z = [1, Z 2 , ■ ■ ■, Z 2 m- P -i] a non-zero vector of Lagrange multipliers (we let Z\ — 1 
w.l.o.g.). One obtains that i is constituted by 

• 2m — p — 1 polynomials of multidegree bounded by (1,1, 0) with respect to (x, y, z), 
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• n — 1 polynomials of multidegree bounded by (0, 1 , 1) with respect to (x, y, z), 

• p polynomials of multidegree bounded by (1, 0,1) with respect to (x, y, z), 

that is by n + 2m — 2 polynomials in n + 2 m — 2 variables. 


Lemma 5 With the above notations, the number of isolated solutions of Z(£) is at most 

'2m — p — l\ / n — 1 \ (p\ 


S(m,n,p ) = 


7i — £ I \2m — 2p — 2 + 


where £ 6 {max{0, n — 2m + p + 1},..., min{p, n — 2m + 2p + 1}}. 


Proof : By [20, Proposition 11.1], this degree is bounded by the multilinear Bezout bound 
S(m,n,p ) which is the sum of the coefficients of 

(s x + Sy) 2m P I (s y + S^) n 1 (s x + S Z ) P G Q[s x , Sy, S 2 ] 

modulo / = ^s" +1 , s p+1 , The conclusion comes straightforwardly by technical 

computations. □ 

With input £, the homotopy system t is constituted by 2m—p — l, 7 i—l and p polynomials 
of multidegree respectively bounded by (1,1, 0,1), (0,1,1,1) and (1, 0,1,1) with respect 
to (x, y, z ,t). We prove the following. 

Lemma 6 degZ(t) e 0(pn(2m — p)S(m,n,p)). 

Proof of Lemma 6: We use Multilinear Bezout bounds as in the proof of Lemma 5. 
The degree of Z(t) is bounded by the sum of the coefficients of 

(s x + s y + s t ) 2m p 1 (s y + s z + s t ) n 1 (s x + s z + s t ) p 

modulo I = <^s" +1 , s p+1 , C Q[s x , s y , s z , s*]. Since the variable s t can appear 

up to power 1, the previous polynomial is congruent to P\ + P 2 + P 3 + P 4 modulo /, with 

Pi = (s x + Sy) 2m p 1 (Sj / + S^) n 1 (s a; + Sj / ) p 

P 2 = (2m -p - l)si(s x + s y ) 2m - p ~ 2 (s y + S 2 ) n_1 (s x + s z ) p , 

P 3 = (n - 1 )s t (s y + s z ) n ~ 2 (s x + s y ) 2m_p_1 (s 3; + s.) p 
P 4 = P St(s x + s z ) p 1 (s x + s y ) 2m p 1 (s y + s^) n b 

We denote by A (Pf) the contribution of P* to the previous sum. 

Firstly, observe that A(P x ) = 5(m,n,p) (compare with the proof of Lemma 5). Defining 
Xi = max{0, n — 2m + p + 1} and y 2 = min {p,n — 2m + 2p + 1}, one has A (Pi) = 

S(m,n,p) = E£ n 7(<) with 7 (<) = (^T') LV- 2 +*) (?)■ 
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Write now P 2 = (2 m — p — 1 )s t P 2 , with P 2 G Q[a;, y,z]. Let A(P 2 ) be the contribution 
of P 2 , that is the sum of the coefficients of P 2 modulo I' = ^s" +1 , s^ +1 , so that 

A(P 2 ) = (2 m — p — 1)A(P 2 ). Then 


A(A) = £ 

i,j/ 


f 2m — p — 2 n 


n 



where the sum runs in the set defined by the inequalities 


i + £ < n, 2 m — p — 2 — i + j < p, n — 1—j+p — £< 2m — p — 2. 


Now, since P 2 is homogeneous of degree n + 2m — 3, only three possible cases hold: 

Case (A), i + £ = n, 2m — p — 2 — i + j = p and n — 1— j +p — £ = 2m — p — 3. Here the 
contribution is S a = J2el ai Pa(£) with 


^o(f) 


/ 2m — p — 2\ / n — 1 \ £p\ 

V n-i J {2m-2p~3 + £j {i) 1 


and = max{0, n — 2m + p + 2}, a 2 = min{p, n — 2m + 2p + 2}. Suppose first that £ is 
an admissible index for A(Px) and S a , that is max-fy!, ax} = ay < £ < Xi — m in{x 2 ,ct 2 }- 
Then: 


2m — p — 1 

n — £ 


<Pa{£) < 

= (*) 7 (*) 


n — 1 \ /p\ 

v 2m — 2p — 3 + “ 

2m — 2p — 2 + 


with T(£) = 


n — (2 m — 2p — 2 + 


The rational function £ i —> V L(£) is piece-wise monotone (its first derivative is positive), 
and its unique possible pole is £ — n — 2m + 2p + 2. Suppose that this value is a pole 
for 'L(f'). This would imply a 2 = n — 2m + 2p + 2 and so y 2 = n ~ 2m + 2p + 1; since 
£ is admissible for A(Px), then one would conclude a contradiction. Hence the rational 
function T(£) has no poles, its maximum is atteined in y 2 and its value is T(y 2 ) = n — 1. 
Hence (p a (£) A (n — 1)7(^). Now, we analyse any possible case: 


(Al) yx = 0, aq = 0. This implies y 2 = n — 2m + 2p + 1, a 2 — n — 2m + 2p + 2. We 
deduce that 


X2 X2 

Va{£) + (fa(a 2 ) < (n - 1) £ j(£) + 

e=o e=o 

+ Pa(a 2 ) < (n - l)A(Pi) + (p a (a 2 ). 

In this case we deduce the bound 5 a < nA(Px). 

(A2) yx = 0, ai = n — 2m + p + 2. This implies X 2 = n — 2m + 2p + 1, a 2 = p. In this 
case all indices are admissible, and hence we deduce the bound 5 a < (n — l)A(Px). 

(A3) xi — n ~ 2m + p + 1. This implies a\ — n — 2m + p + 2, X 2 = P, ol 2 = p. Also in 
this case all indices are admissible, and 8 a < (n — l)A(Px). 
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Case (B). i + £ — n, 2m — p — 2 — i + j = p — 1 and n — 1 — j + p — £ = 2m — p — 2. Here 
the contribution is 5b = Tb(£) where 


MO 


(2m — p — 2\ / n — 1 \ fp\ 

V n - £ J\2m-2p-2 + e)[e)' 


One gets 5b < A (Pi) since the sum above is defined in max{0, n2m + p + 2} < £ < 
min{p, n — 2m + 2p + 1}, and the inequality (fb(0 < 7(^) holds term-wise. 

Case (C) i + £ = n — 1, 2m — p — 2 — i+j = p and n— 1— j + p — £ — 2m — p — 2. Here 
the contribution is 5 C = ^c(^) where 


MO 


( 2m — p — 2\ / n — 1 \ /p\ 

\n-\ -£ J \2m -2p-2 + £J\i)' 


One gets 5 C < A (Pi) since the sum above is defined in max{0, n2m + p + 1} < £ < 
min {p,n — 2m + 2 p + 1}, and the inequality ip c (£) < q(f') holds term-wise. We conclude 
that 5 a < nA(Pi), 5b < A(Pi) and 5 C < A (Pi). Hence A (P 2 ) = (2m—p — l)(5 a + 5b + 5 c ) G 
0(n(2m — p)A(Pi)). 

Analogously to A(P 2 ), one can conclude that A(P 3 ) e 0(n(n + 2?n — p)A(P 1 )) and 
A(P 4 ) G 0(pn(n + 2m — p)A(Pi)). □ 

Estimates. 

We provide the whole complexity of ZeroDimSolveMaxRank. 


Theorem 7 Let 5 = 5(m,n,p ) 6e given by Lemma 5. Then ZeroDimSolveMaxRank with 
input £(H m , Up, Vp) computes a rational parametrization within 

0~(p{n + 2m)° ( A(2m — p)5' 2 ), 


arithmetic operations over Q. 

Proof : The polynomial entries of the system t (as defined in the previous section) are 
cubic polynomials in n + 2m — 1 variables, so the cost of their evaluation is in 0((n + 
2m) 3 ). Applying [15, Theorem 5.2] and bounds given in Lemma 5 and 6 yield the claimed 
complexity estimate. □ 

From Lemma 5, one deduces that for all 0 < p < r, the maximum number of complex 
solutions computed by ZeroDimSolveMaxRank is bounded above by 5(m,n,p). We deduce 
the following result. 


Proposition 8 Let H be a m x m, n — variate linear Hankel matrix, and let r < m — 1. 
The maximum number of complex solutions computed by LowRankHankel with input ( H,r ) 
is 


2 m 


n r 

k=2m—2r p =0 


where 5(m, k,p ) is the bound defined in Lemma 5. 
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Proof : The maximum number of complex solutions computed by ZeroDimSolve is the 
degree of Using, the multilinear Bezout bounds, this is bounded by the 

coefficient of the monomial s™Sy in the expression ( s x J rS y ) 2rn ~ r ~ 1 , that is exactly ( 2m ~ r_1 V 
The proof is now straightforward, since ZeroDimSolveMaxRank runs r + 1 times at each 
recursive step of LowRankHankel, and since the number of variables decreases from n to 
2 m — 2 r. □ 


5 Proof of Proposition 2 

We start with a local description of the algebraic sets defined by our Lagrange systems. 
This is obtained from a local description of the system defining J(H , u, p). Without loss of 
generality, we can assume that u = (0,..., 0,1) in the whole section: such a situation can 
be retrieved from a linear change of the y-variables that leaves invariant the x-variables. 


5.1 Local equations 


Let (x, y) G JC(H,u,p). Then, by definition, there exists a p x p minor of H(x ) that is 
non-zero. Without loss of generality, we assume that this minor is the determinant of the 
upper left p x p submatrix of H. Hence, consider the following block partition 


H(x) 


N Q 
P R 


( 1 ) 


with N G Q[xf x P, and Q G Q[x] p , P G Qfxp™- 2 ?" 1 )^ and R G Q[x] 2m - 2 P-'. We 
are going to exhibit suitable local descriptions of )C(H,u,p) over the Zariski open set 
On C C n+P+1 defined by detiV ^ 0; we denote by Q[x, yjdetiv the local ring of Q[x, y] 
localized by det N. 


Lemma 9 Let N, Q, P, R be as above, and u G Q p+1 — {0}. Then there exist {qi}i<i< p C 
Q[x]det n an d {qi}i<i< 2 m- 2 P -i C Q[x] det 7 v such that the constructible set JC(H,u,p) fl On 
is defined by the equations 

Vi ~ qi(x) =0 1 < i < P 

qfix) =0 1 < i < 2 m — 2p — 1 

V P +1 -i = o. 


Proof : Let c = 2 m — 2p—l. The proof follows by the equivalence 


N Q 
P R 


y = 0 iff 


I P 0 

-P Ic 


N- 1 0 

0 Ic 


N Q 
P R 


y = o 


in the local ring Q[x,y]d e tiv, that is if and only if 


Ip N~ l Q 
0 R - PN~ l Q 


y = 0 
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Recall that we have assumed that u = (0,..., 0,1); then the equation uy = 1 is y p+ 1 = 1. 
Denoting by g* and g* respectively the entries of vectors —N~ X Q and — (R — PN~ l Q) ends 
the proof. □ 

The above local system is denoted by f G Q[x, y]^jv- Th e Jacobian matrix of this 
polynomial system is 

Dt= \ ” '' ,° 

with q = (gi(x),..., g 2m _ 2 P -i(x)). Its kernel defines the tangent space to /C(L7, u,p) HOn- 
Let w = (wi,..., w„) G C" be a row vector; we denote by 7r w the projection 7r w (x, y) = 
wiXi + • • • + w n x n . Given a row vector v G C 2m_p+1 , we denote by wlagrange(f, v) the 
following polynomial system 


fj (§) h) [^lj ■ ■ ■ j %2rri—pi ^2m— p+l] 


D f 
w 0 


vz 


( 2 ) 


For all 0 < p < r, this polynomial system contains n + 2m + 2 polynomials and n + 2m + 2 
variables. We denote by L p (f, v, w) the set of its solutions whose projection on the (x, y)- 
space lies in On- 

Finally, we denote by wlagrange(f, v) the polynomial system obtained when replacing f 
above with f = f(L7, u, p). Similarly, its solution set is denoted by L p (f, v, w). 


5.2 Intermediate result 

Lemma 10 Let YZ C £,( 2m - r )( n + 1 ) the non-empty Zariski open set defined by Proposi¬ 
tion 1, H e and 0 < p < r. There exist non-empty Zariski open sets Y C C 2m_p and 
f cC" such that if v e Y and w e W, the following holds: 

(a) the set C p ( f,v,w) = £(f, v, w) n {(x,y, z) | ranklL(x) = p} is finite and the 
Jacobian matrix of wlagrange(f, v) has maximal rank at any point of C p ( f, v, w); 

(b) the projection of C p { f, v, w) in the (x, y) -space contains the critical points of the 
restriction of 7r w restricted to K,(H,u,p). 

Proof : We start with Assertion (a). 

The statement to prove holds over IC(H, u,p); hence it is enough to prove it on any open 
set at which one p x p minor of H is non-zero. Hence, we assume that the determinant of 
the upper left p x p submatrix of H is non-zero; On C C n+P+1 is the open set defined 
by det N ^ 0, and we reuse the notation introduced in this section. We prove that there 
exist non-empty Zariski open sets Yfj C C 2m ~ p and Wn C C n such that for v G Yjf and 
w G Wn-, £ p (f, v, w) is finite and that the Jacobian matrix associated to wlagrange(f, v) 
has maximal rank at any point of £ p (f, v, w). The Lemma follows straightforwardly by 
defining Y' (resp. W) as the intersection of Yfj (resp. Wn) where N varies in the set of 
p x p minors of Hfx). 

Equations h yield Zj = 0 for j = 2 m — 2 p, ..., 2 m — p, and can be eliminated together 
with their z variables from the Lagrange system wlagrange(f, v). It remains z-variables 
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Zi,, z 2m - 2 p-i, Zim-p+i'i we denote by fl C C 2m 2p the Zariski open set where they don’t 
vanish simultaneously. 

Now, consider the map 


q : OjvxOx C n —X C n+2m - p 

( x , y, z, w) i—x (f,g) 

and, for w G C n , its section map g w (x, y, z) = g(x, y, z,w). We consider v G C 2m_p and 
we denote by z the remaining z— variables, as above. Hence we define 

Q : 0 N XflxCx C 2m " 2p -» C n+ 2 m-p+l 

(x,y,z,w,v) i—X (f,g,v'z-l) 

and its section map Q w ,v(x, y,z) = g(x, y, z, w, v). We claim that 0 G C n+2m ~ p (resp. 
0 G C n+2m ~ p+1 ) is a regular value for q (resp. Q). Hence we deduce, by Thom’s Weak 
Transversality Theorem, that there exist non-empty Zariski open sets Wn C C n and 
Y N c C 2m ~ 2p such that if w G Wn and v G Yn, then 0 is a regular value for q w and Q w ,v- 

We prove now this claim. Recall that since H G M 1 , the Jacobian matrix D x y f has 
maximal rank at any point (x, y) G Z(i\ Let (x, y, z, w) G (^(O) (resp. (x, y, z, w, v) G 
Q _1 (0)). Hence (x, y) G Z(f). We isolate the square submatrix of D g(x, y, z, w) obtained 
by selecting all its rows and 

• the columns corresponding to derivatives of x, y yielding a non-singular submatrix 
of D X y f(x, y); 

• the columns corresponding to the derivatives w.r.t. wy,..., w n , hence this yields a 
block of zeros when applied to the lines corresponding to f and the block I n when 
applied to g. 

For the map Q, we consider the same blocks as above. Moreover, since (x, y, z,w, v) G 
Q _1 (0) verifies v'z — 1 = 0, there exists ^ such that zg ^ 0. Hence, we add the derivative 
of the polynomial v'z — 1 w.r.t. vg, which is zg ^ 0. The claim is proved. 

Note that q , vv 1 (0) is defined by n+2m—p polynomials involving n+2m— p+1 variables. We 
deduce that for w G #at, q'“ 1 (0) is either empty or it is equidimensional and has dimension 
1. Using the homogeneity in the z-variables and the Theorem on the Dimension of Fibers 
[21, Sect. 6.3, Theorem 7], we deduce that the projection on the (x, y)-space of (^(O) 
has dimension < 0. We also deduce that for w G Wn and v G Yn, Qw'v(O) is either empty 
or finite. 

Hence, the points of Qvw(O) are in bijection with those in £(f, v, w) forgetting their 
0-coordinates corresponding to Zj = 0. We define = Yn xPc C 2m ^ 2p . We deduce 
straightforwardly that for v G Y^ and w G the Jacobian matrix of wlagrange(f, v) 
has maximal rank at any point of C p ( f, v, w). By the Jacobian criterion, this also implies 
that the set £ p (f, v, w) is finite as requested. 

We prove now Assertion (b). 
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Let f cC" and Y' C C 2m_p be the non-empty Zariski open sets defined in the proof of 
Assertion (a). For w G W and v G Y' , the projection of £ p (f, v, w) on the (x, y)— space 
is finite. Since H G JC(H,u,p) is smooth and equidimensional. 

Since we work on IC(H, u, p), one of the pxp minors of iL(x) is non-zero. Hence, suppose 
to work in OnOJC(H, u ,p) where On C C n+P+1 has been defined in the proof of Assertion 
(a). Remark that 

crit ( 7 r w , 1C(H, u,p)) = IJ crit ( 7 r w , O n fl JC(H, u,p)) 

N 

where N runs over the set of pxp minors of iL(x). We prove below that there exists a 
non-empty Zariski open set Y C C 2m_p such that if v G Y, for all N and for w G W, the 
set crit ( 7 r w , On H JC(H, u,p)) is finite and contained in the projection of C p ( f, v, w). This 
straightforwardly implies that the same holds for crit ( 7 r w , /C(L/j u,p)). 

Suppose w.l.o.g. that N is the upper left pxp minor of iL(x). We use the notation f, g, h 
as above. Hence, the set crit ( 7 t w ,Ojv H fC(H,u,p)) is the image by the projection 7r xy 
over the (x, y)— space, of the constructible set defined by f, g, h and z ^0. We previously 
proved that, if w G Wn, <? _1 (0) is either empty or equidimensional of dimension 1. Hence, 
the constructible set defined by f, g, h and z 7 ^ 0, which is isomorphic to g _ 1 ( 0 ), is either 
empty or equidimensional of dimension 1 . 

Moreover, for any (x, y) G crit (7T W , On H /C(//, u,p)), 7 r“y(x, y) has dimension 1, by the 
homogeneity of polynomials w.r.t. variables z. By the Theorem on the Dimension of 
Fibers [21, Sect. 6.3, Theorem 7], we deduce that crit ( 7 r w , On H IC(H, u,p)) is finite. 

For (x, y) G crit ( 7 r w , OnO)C(H, u ,p)), let ^ x ,y ),n C C 2m_p be the non-empty Zariski open 
set such that if v G 'f'^ y ) t N the hyperplane v'z — 1 = 0 intersects transversely 7 r“ y (x, y). 
Recall that C C 2m-p has been defined in the proof of Assertion (a). Define 

Vn = Vn n fl ),N 
(x.y) 

and V = Hat ^n- This concludes the proof, since Y is a finite intersection of non-empty 
Zariski open sets. □ 


5.3 Conclusion 


We denote by ,M\ C GL(n, C) the set of non-singular matrices M such that the first row 
w of M _1 lies in the set W given in Lemma 10: this set is non-empty and Zariski open 
since the entries of M~ l are rational functions of the entries of M. Let Y C C 2m_p be 
the non-empty Zariski open set given by Lemma 10 and let v G Y. Let ei be the row 
vector (1, 0,..., 0) G Q n and for all M G GL(n, C), let 


M 


M 0 
0 I m 


Remark that for any M G ^#1 the following identity holds: 


' Df(H M ,u,p) ' 


' Df(H,u,p) ' 

ei 0 ••• 0 


w 0 • • • 0 


19 








We conclude that the set of solutions of the system 


f z 


Df(H,n,p) 

w 0 • • • 0 


\'z — 1 


(3) 


is the image by the map (x, y) H > M 1 (x, y) of the set S of solutions of the system 




Df(H,n,p) 

ei 0 ••• 0 


v'z — 1 


(4) 


Now, let <p be the projection that eliminates the last coordinate Z 2 m -p+i ■ Remark that 
<P{S) = L p (f M , v,ei). 

Now, applying Lemma 10 ends the proof. 

□ 


6 Proof of Proposition 3 

The proof of Proposition 3 relies on results of [14, Section 5] and of [19]. We use the same 
notation as in [14, Section 5], and we recall them below. 


Notations For Z C €" of dimension d, we denote by fli(Z) its i—equi dimensional 
component, i — 0,.. ., d. We denote by S^{Z) the union of: 


• VLq(Z) U • • • U PL d -i(Z) 

• the set sing (f ld(Z)) of singular points of fl d (Z). 

Let 7q be the map (xi,..., x n ) —>■ (xi,..., Xj). We denote by ^( 7 q, Z) the Zariski closure 
of the union of the following sets: 


• £Iq{Z) U • • ■ U Qj_i(Z); 

• the union for r > i of the sets crit ( 7 q, reg {VL r (Z))). 

For M G GL(n, C) and Z as above, we define the collection of algebraic sets {Oi(Z M )} 0 <i< d 
as follows: 

• O d (Z M ) = z M - 

. O^Z™) = y{o l+l {z M )) u ^(tt j+1 , o l+ 1 (z M )) u 

U e jf(7Ti + i, Z AI ) for i — 0,..., d — 1. 

We finally recall the two following properties: 

Property P (Z). Let 2 C €" be an algebraic set of dimension d. We say that M e GL(n, C) 
satisfies P (Z) when for alii = 0,1,..., d: 
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1. Oi{Z AI ) has dimension < i and 

2. Oi{Z AI ) is in Noether position with respect to xi,..., X{. 

Property Q. We say that an algebraic set Z of dimension d satisfies Qi(Z) (for a given 
1 < i < d) if for any connected component C of Z flM" the boundary of 7 Tj(C) is contained 
in 7Ti(Oi-i(Z) n C). We say that Z satisfies Q if it satisfies Qi,..., Q d- 

Let Z C C" be an algebraic set of dimension d. By [14, Proposition 15], there exists a 
non-empty Zariski open set C GL(n, C) such that for M G D GL(n, Q) Property 
P (Z) holds. Moreover, if M G GL(n, Q) satisfies P (Z), then Qi{Z M ) holds for i — 1,..., d 
[14, Proposition 16]. We use these results in the following proof of Proposition 3. 

Proof : We start with assertion (a). Let ^#2 C GL(n, C) be the non-empty Zariski open 
set of [14, Proposition 17] for Z = Pip. for M G ^ 2 , M satisfies P (PL p ). Remark that the 
connected components of Pi p D R n and are in bijection with those of Pi AI fl M n (given by 
C GG C M ). Let C M C Pip 1 fl M n be a connected component of PL^ fl M n . Let 7Ti be the 
projection on the first variable 7Ti: M n —> M, and consider its restriction to 'H( v/ nR n . Since 
M G ^# 2 , by [14, Proposition 16] the boundary of ni(C M ) is included in TTi( y O 0 (Pi^f) C\C M ) 
and in particular in ni(C M ). Hence tti(C ai ) is closed. 

We prove now assertion (b). Let M G ^# 2 , C a connected component of PL p fl M" and 
aGlbe in the boundary of 7 Ti(C' m ). By [14, Lemma 19] 7r j" 1 (a) fl C M is finite. 

We claim that, up to genericity assumptions on u G Q p+1 , for x G 7rj“ 1 (tt) fl C M , the 
linear system y H- i(H M ,u,p) has at least one solution. We deduce that there exists 
a non-empty Zariski open set '^c , ,x C C p+1 such that if u G ^c,x H Q p+1 , there exists 
y G Q p+1 such that (x, y) G J(H M , u,p). One concludes by taking 

^ = n n ^x, 

CcW P nr x ^-i( Q) nc« 

which is non-empty and Zariski open since: 

• the collection {C C PL p fl M n connected component} is finite; 

• the set 7 rj _1 (a) fl C M is finite. 

It remains to prove the claim we made. For x G 7rf 1 (a) fl C M , the matrix LT(x) is rank 
defective, and let p' < p be its rank. The linear system 

' H(x) 

u 

has a solution if and only if 


rank 

' H(x) ‘ 

= rank 

' H (x) 0 " 


u 


u 1 


and the rank of the second matrix is p' + 1. Denoting by C C p+1 the complement in 
C p+1 of the p '—dimensional linear space spanned by the rows of if(x), proves the claim 
and concludes the proof. □ 


0 

1 
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7 Experiments 


The algorithm LowRankHankel has been implemented under Maple. We use the FGb 
[ 9] library implemented by J.-C. Faugere for solving solving zero-dimensional polynomial 
systems using Grobner bases. In particular, we used the new implementation of [10] for 
computing rational parametrizations. Our implementation checks the genericity assump¬ 
tions on the input. 

We test the algorithm with input m x m linear Hankel matrices R(x) = Hq + X\H\ + 
... + x n H n , where the entries of Ho,, H n are random rational numbers, and an integer 
0 < r < m — 1. None of the implementations of Cylindrical Algebraic Decomposition 
solved our examples involving more that 3 variables. Also, on all our examples, we found 
that the Lagrange systems define finite algebraic sets. 

We compare the practical behavior of LowRankHankel with the performance of the library 
RAGlib, implemented by the third author (see [17]). Its function PointsPerComponents, 
with input the list of (r +1)—minors of R(x), returns one point per connected component 
of the real counterpart of the algebraic set H r , that is it solves the problem presented in 
this paper. It also uses critical point methods. The symbol oo means that no result has 
been obtained after 24 hours. The symbol matbig means that the standard limitation in 
FGb to the size of matrices for Grobner bases computations has been reached. 

We report on timings (given in seconds) of the two implementations in the next table. 
The column New corresponds to timings of LowRankHankel. Both computations have 
been done on an Intel(R) Xeon(R) CPU F7540 @2.00GHz 256 Gb of RAM. We remark 
that RAGlib is competetitive for problems of small size (e.g. m = 3) but when the 
size increases LowRankHankel performs much better, especially when the determinantal 
variety has not co-dimension 1. It can tackle problems that are out reach of RAGlib. 
Note that for fixed r, the algorithm seems to have a behaviour that is polynomial in nm 
(this is particularly visible when m is fixed, e.g. to 5). 

Finally, we report in column Total Deg the degree of the rational parametrization obtained 
as output of the algorithm, that is the number of its complex solutions. We observe that 
this value is definitely constant when m, r are fixed and n grows, as for the maximum 
degree (column MaxDeg) appearing during the recursive calls. 

The same holds for the multilinear bound given in Section 4 for the total number of 
complex solutions. 
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(m, r, n) 

RAGlib 

New 

TotalDeg 

MaxDeg 

(3,2,2) 

0.3 

5 

9 

6 

(3,2,3) 

0.6 

10 

21 

12 

(3,2,4) 

2 

13 

33 

12 

(3,2,5) 

7 

20 

39 

12 

(3,2,6) 

13 

21 

39 

12 

(3,2,7) 

20 

21 

39 

12 

(3,2,8) 

53 

21 

39 

12 

(4,2,3) 

2 

2.5 

10 

10 

(4,2,4) 

43 

6.5 

40 

30 

(4,2,5) 

56575 

18 

88 

48 

(4,2,6) 

(X) 

35 

128 

48 

(4,2,7) 

(X) 

46 

143 

48 

(4,2,8) 

oo 

74 

143 

48 

(4,3,2) 

0.3 

8 

16 

12 

(4,3,3) 

3 

11 

36 

52 

(4,3,4) 

54 

31 

120 

68 

(4,3,5) 

341 

112 

204 

84 

(4,3,6) 

480 

215 

264 

84 

(4,3,7) 

528 

324 

264 

84 

(4,3,8) 

2638 

375 

264 

84 

(5,2,5) 

25 

4 

21 

21 

(5,2,6) 

31176 

21 

91 

70 

(5,2,7) 

x 

135 

199 

108 

(5,2,8) 

oo 

642 

283 

108 

(5,2,9) 

oo 

950 

311 

108 

(5,2,10) 

oo 

1106 

311 

108 

mmm 

2 

2 

20 

20 

m 

202 

18 

110 

90 

HI § 

(X 

583 

338 

228 

IgBBl 

(X 

6544 

698 

360 

Hi 

X 

28081 

1058 

360 

ran 

X 

X 

- 

- 


Tabic 1: Timings and degrees 
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(m, r, n ) 

RAGlib 

New 

TotalDeg 

MaxDeg 

HSE1 

1 

5 

25 

20 

rail 

48 

30 

105 

80 


8713 

885 

325 

220 

m 

oo 

15537 

755 

430 

RH 

oo 

77962 

1335 

580 

(6,2,7) 

oo 

6 

36 

36 

(6,2,8) 

oo 

matbig 

- 

- 

(6,3,5) 

oo 

10 

56 

56 

(6,3,6) 

oo 

809 

336 

280 

(6,3,7) 

oo 

49684 

1032 

696 

(6,3,8) 

oo 

matbig 

- 

- 

(6,4,3) 

3 

5 

35 

35 

(6,4,4) 

oo 

269 

245 

210 

(6,4,5) 

oo 

30660 

973 

728 

(6,4,6) 

oo 

oo 

- 

- 

(6,5,2) 

1 

9 

36 

30 

(6,5,3) 

915 

356 

186 

150 

(6,5,4) 

oo 

20310 

726 

540 

(6,5,5) 

oo 

oo 

- 

- 


Table 2: Timings and degrees (continued) 
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